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Abstract. We study the single spin asymmetries for the D-meson production, Drell-Yan lepton-pair 
production and the direct-photon production in the pp colHsion induced by the twist-3 three-gluon 
correlation functions in the transversely polarized nucleon. We present a corresponding polarized 
cross section formula in the leading-order with respect to the QCD coupling and a model calculation 
for the asymmetries, illustrating the sensitivity to the form of the correlation functions. 
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1. INTRODUCTION 

Understanding the origin of the large single spin asymmetries (SSAs) observed in 
various high-energy semi-inclusive processes have been a big challenge during the past 
decades. The SSA can be generated as a consequence of the multiparton correlations 
inside the hadrons in the coUinear factorization approach which is valid when the 
transverse momentum of the particle in the final state can be regarded as hard. Among 
such multiparton correlations, purely gluonic correlations have a potential importance, 
since gluons are ample in the nucleon. The best way to probe such gluonic correlations is 
the measurement of SSA for a heavy meson production in semi-inclusive deep inelastic 
scattering (SIDIS) and the pp collision [1, 2], since the heavy-quark pair one of which 
fragments into the final-state meson is mainly produced, respectively, by the photon- 
gluon and gluon-gluon fusion mechanisms. The measurement of SSA for the D-meson 
production in the pp collision is ongoing at RHIC [3]. 

In this report, we study the contribution of the three-gluon correlation functions 
representing such multigluonic effects to SSA in the pp collision for the D-meson 
production (p^p — DX) [4], Drell-Yan lepton-pair production (p^p — ?■ 7*X) and the 
direct-photon production (p^p — t- 7X')[5]. We will present the corresponding single- 
spin dependent cross sections by applying the formalism developed for SIDIS, ep^ — t- 
eDX [2]. We will also present a model calculation of the SSAs for p^p DX and 
p^p — >■ induced by the three-gluon correlation functions in comparison with the 
RHIC preliminary data for the former [3]. 



2. THREE-GLUON CORRELATION FUNCTIONS 



Three-gluon correlation functions in the transversely polarized nucleon are defined as 
the color-singlet nucleon matrix element composed of the three gluon's field strength 
tensors F"^. Corresponding to the two structure constants for the color SU(3) group, 
di,ca and /^^a. one obtains two independent three-gluon correlation functions 0{x\,X2) 
and A^(xi,X2) as [2] 
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where S is the transverse-spin vector for the nucleon, n is the light-like vector satisfying 
p n = 1, and we have used the shorthand notation as F^" = F^^np etc. The gauge-link 
operators which restore gauge invariance of the correlation functions are suppressed in 
(1) and (2) for simplicity. The nucleon mass is introduced to define 0{xi,X2) and 
N{xi,X2) dimensionless . 



3. D-MESON PRODUCTION IN COLLISION 

Applying the formalism for the contribution of the three-gluon correlation functions to 
SSA developed in [2], the twist-3 cross section for p^{p,S±^) + p{p') — )■ D{Ph) +X with 
the center-of-mass energy can be obtained in the following form [4]: 
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where 5c = I and 5c = —1, Df{z) represents the c — > D or c — > D fragmentation 

functions, G{x') is the unpolarized gluon density, pc is the four-momentum of the c 
(or c) quark (mass nic) fragmenting into the final D (or D) meson and s, t, u are 
defined as 5 = {xp -\- x' p')^ , t — (xp — pc)^ — m^, u = {x'p' — pcY — m^. The hard cross 
sections o-<5i,02,03,04 ^m,N2,m,N4 listed in [4]. The cross section (3) receives 
contributions from 0{x,x), O{x,0), N{x,x) and A^(x,0) separately, which differs from 
the previous result [1]. 

We perform numerical estimate for A/^ based on (3). For the RHIC kinematics, 
we found that the terms with ^03,04,N3,N4 negligible compared with those with 
^oi,02,Ni,N2 ^^Yat o^^ ^ d*^^ ~ d^^ ^ —a^^. One can thus regard the cross section 
as a function of the correlation functions 0{x) = 0{x,x) + 0{x, 0) and N{x) = N{x,x) — 



N{x,0) to a very good approximation. We assume the relation 0{x) = ±N{x) together 
with 0{x,x) = O{x,0) and N{x,x) = —N{x,0) for simplicity. For the functional form of 
each function, we employ the following two models: 

Model 1 : 0{x) = KgxG{x), (4) 
Model 2 : 0{x) = isT^ v^G(jc) , (5) 

where Kq and K'q are the constants to be determined so that the calculated asymmetry is 
consistent with the RHIC data [3]. 

For the numerical calculation, we use GJR08 [6] for G{x) and KKKS08 [7] for 
We calculate Ajv for the D and D mesons at the RHIC energy at \fS = 200 GeV and the 
transverse momentum of the D-meson Pj = 2 GeV. We set the scale of all the distribution 

and fragmentation functions at jU = w + with the charm quark mass = 1 .3 GeV. 




FIGURE 1. (a) for and (b) An for for Model 1 in (4) with 0{x) ^ N(x) and Kq = 0.004, and 
(c) for D° and (c) An for D° for Model 2 in (5) with 0{x) = N{x) and /T^ = 0.001. Short bars denote 
the RHIC preliminary data taken from [3]. 

Fig. 1 shows the result of A^y for the D° and D*^ mesons with the relation 0{x) = 
N{x) together with the preliminary data [3] denoted by the short bars. The sign of 
the contribution from 0{x) changes between and as shown in (3) and works 
constructively (destructively) for D*^ 0^) for the case 0{x) — N{x). The values Kq = 
0.004 and Kq = 0.001 have been determined so that Apj does not overshoot the RHIC 
data. If one adopts the relation 0{x) = —N{x), the result for the and mesons 
will be interchanged. The rising behavior of Aj^ at large xp as shown in Figs. 1(a) 
and (c) is originated from the derivative of 0{x) and N{x), as in the case of the soft- 
gluon-pole (SGP) contribution for the quark-gluon correlation function. By comparing 
the results for the models 1 and 2 in Figs. 1(a) and (c), one sees that the behavior of 
AjM at Xf < depends strongly on the small-x behavior of the three-gluon correlation 
functions. Therefore Aa? at x/r < is useful to get constraint on the small-jc behavior of 
the three-gluon correlation functions. 

4. DRELL-YAN LEPTON-PAIR PRODUCTION IN pp COLLISION 

The analysis of the previous section can be extended to A^ for the Drell-Yan lepton-pair 
production induced by the three-gluon correlation function. The corresponding twist- 
3 diagrams for the hard part are shown in Fig. 2, which give rise to SGP contribution 
at xi = X2 due to the initial-state interaction between the extra incoherent gluon from 



the polarized proton and the quark coming out of the unpolarized proton. From these 




FIGURE 2. The diagrams for the three-gluon contribution to the twist-3 cross section for p^p 'fX. In 
the coUinear Umit, the momenta k\ and coming from the polarized nucleon are set to ki ~ xip (/ = 1 , 2) 
and the pole contribution from the bared propagator gives rise to the SSA. 

diagrams, one obtains for the single-spin-dependent cross section for the Drell-Yan 
process, p^{^p, 5^) + — > 7*(^) with the invariant mass = for the lepton- 
pair as [5] 
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where y is the rapidity of the virtual photon, Uem — 1/137 is the QED coupling constant, 
and the partonic hard cross sections are defined as 
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with the number of colors N = 3 and s = {xp -\-x' p'Y, t = [xp — <?)^ and u = [x'p' — q)^. 
For a large Q^, di differs from 02 significantly, and 4 are not negligible. Therefore the 
cross section (6), in general, depends on the four functions 0{x,x), O{x,0), N{x,x) and 



N{x, 0) independently unlike the case for p^p — )■ DX in the previous section, where the 
twist-3 cross section can be regarded as a function of the combination 0(jc, jc) + 0(jc, 0) 
andA^(;c,;c) -A^(;c,0). 

The An for the Drell-Yan process receives contribution not only from the three-gluon 
correlation functions but also from the quark-gluon correlation functions, for which the 
twist-3 cross section have been derived in [8, 9, 10]. The sum of (6) and those from the 
quark-gluon correlation functions gives the complete leading-order cross section for the 
asymmetry. 



5. DIRECT PHOTON PRODUCTION IN pp COLLISION 

Taking the ^ — )■ limit of (6), the twist-3 cross section for the direct photon production, 
P^{P-iS±) + p{p') — )■ y{q) +X, induced by the three-gluon correlation functions can be 
obtained as [5] 
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where fa{^) is the twist-2 unpolarized quark density and 5a = 1(— 1) for quark (an- 
tiquark). As shown in (8), the combinations 0{x) = 0{x,x) + O{x,0) and N{x) = 
N{x,x) —N{x,0) appear in the cross section accompanying the common partonic hard 
cross section which is the same as the twist-2 hard cross section for the qg qy scat- 
tering. The origin of this simplification can be clearly understood in terms of the master 
formula for the contribution from the three-gluon correlation functions developed in 
[1 1, 4, 5]. We also note that the above result (8) differs from the previous study in [12]. 

We have performed a numerical calculation for A]^ based on the models used for 
p^p — )■ DX in Sec. 3. For the models 1 and 2, we calculate the asymmetry A]^ for the two 
cases: Case 1 ; 0{x) = A^(jc) and Case 2; 0{x) = -A^(jc). We use GJR08 [6] for fq{x') and 
the models (4) and (5) with Kq = 0.004 and K'^ = 0.001 which are consistent with the 
RHIC data. We calculate A J at the RHIC energy atVS = 200 GeV and the transverse 
momentum of the photon qr = 2 GeV, setting the scale of all the distribution functions 
at II = qj. 

Fig. 3 shows the result for A]^ for each case. One can see Ayy aixp>Q become almost 
zero regardless of the magnitude of the three-gluon correlation functions, while A]^ at 
xp <0 depends strongly on the small-.)c behavior of the three-gluon correlation functions 
as in the case of p^p — > DX. Even though the derivatives of 0{x) and N{x) contribute, Aa? 
is tiny at jcf > due to the small partonic cross section. At xp < 0, large-jc' region of the 
unpolarized quark distributions and the small-x region of the three-gluon distributions 
are relevant. For the above case 1, only antiquarks in the unpolarized nucleon are active 
and thus lead to small aJ^ as shown in Figs. 3(a) and (b). On the other hand, for the 
case 2, quarks in the unpolarized nucleon are active and thus lead to large aJ^ as shown 




FIGURE 3. (a) An for Case 1 with Model 1. (b) An for Case 1 with Model 2. (c) An for Case 2 with 
Model 1. (d) An for Case 2 with Model 2. 



in Figs. 3(c) and (d). Therefore Aj^ at < for the direct photon production could 
provides us with an important information on the relative sign between 0{x) and N{x). 

To summarize, we have studied the SSA for p^p DX, p^p — )• j^X and p^p — )• yX in- 
duced by the three-gluon correlation functions in the polarized nucleon. Combined with 
the known result for the contribution from the quark-gluon correlations, this complete 
the leading-order twist-3 cross sections for these processes. We have also presented a 
model calculation for the asymmetry for the p^ p — )■ DX and p^p — )■ yX at the RHIC en- 
ergy, showing the sensitivity of the asymmetry to the form of the three gluon-correlation 
functions. 
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